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TECHNICAL NOTE 3125 


A SIMPLE MECHANICAL ANALOGUE FOR STUDYING THE DYNAMIC 
STABILITY OF AIRCRAFT HAVING NONLINEAR 
MOMENT CHARACTERISTICS 
By Thomas N. Canning 


SUMMARY 


'Hie analogy between a ball rolling over a contoured surface and a 
pitching and yawing aircraft is developed. Test results from a model 
representing a missile with linear moment characteristics are presented 
to verify the analogy. Several examples of the behavior of nonlinear 
systems are also given. These examples include results from ballistic- 
range firings as' well as analogue tests. 


INTRODUCTION 


Designers of aircraft are frequently confronted with the problem 
of predicting the behavior of aircraft having nonlinear pitching- and 
yawing-moment characteristics. The usual approach, to these problems 
involves tedious calculations or -use of elaborate simulators. A simpler 
approach suitable for some of the problems encountered, particularly 
those of missiles, is the subject of this paper. The approach uses the 
einalogy between a ball rolling in a suitably shaped bowl and a missile 
pitching and yawing in flight. This analogue was devised in order to 
understand a peculiar motion executed by a projectile in the Ames super- 
sonic free-f light range. 


SYMBOLS 


A,B constants defining variation of z with x and y 
a radius of ball, ft 
ag radius of gyration of' ball, ft 

C,D constants defining variation of and Cn with a and p 
Cm 


pitching-moment coefficient of missile 



2 


MCA TN 3125 


Cn yawing-mcanent coefficient of missile 

C(j vector svnn of Cm and Cn for axially symmetric missile 

F force component acting on ball parallel to bowl siirface at point 

of contact which prevents slippage 

g acceleration of gravity, ft/sec^ 

I moment of inertia of ball about a diameter, slug-ft^ 

K constant, ft/sec^ 

Ki, 2 constants, radians/sec^ 

L constant defining rolling friction of ball, l/sec 

I body length, ft 

^1,2 constant defining damping of missile, l/sec 

r y x^ + y^ 

¥ weight of ball, lb 

x,y,z coordinate system for bowl description, x and y horizontal and 
z vertical 

Xjjg/z position of center of gravity as fraction of body length 
a angle of attack, radians 

P angle of sideslip, radians 

6 angle of inclination of bowl surface, radians 

e tan"^ Z! 

X 

0 tan“^ £ 

a 

o V a® + P® 

w rotational velocity of ball, radians/sec 

Superscript 

, differentiation with respect to time 
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BASIC DERIVATICfl? 


Consider the force system acting on a hall rolling on an inclined 
surface . 



X 


Two forces act on the ball to produce translatory acceleration; the 
component of the weight force W parallel to the surface at the point of 
contact, and the tangential force F at the point of contact which pre- 
vents slippage. 


W sin 6 - P = - — ^ 
g cos S 


( 1 ) 


Fa = lu 


Since slippage does not occur 


X = ao) cos 8 
X = ao) cos 6 


Using this relationship 



a cos 6 


Substituting in equation (l) 


sin S 


1 /W I \ X ^ 1 X 
W g cos 6 K cos 8 
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where 


1 

K 




and Bg Is the radius of gyration of the ball. 
If -| 5 [ < 0.15 radian, 


sin 8 cos 8 = sin 8=8 


hz 

Sx 


and it follows that 


X =k5s 
Sx 


The derivation has thus far been limited to motion in the x,z plane. 

In the general case of motion anywhere on the bowl, the component accel- 
erations parallel to the x and y axes, respectively, are: 

i- = K ^ (x,y) 

OX 

y = K (x,y) 

oy 


Similar equations govern the motion of an undamped pitching sind 
yawing missile 

a = Ki Om (a,p) 

P = Ka On (a,p) 


Comparison of these equations shows that a and p are represented 
in the analogy by x and y and that Cm and C^^ are simulated by the 
bowl slopes, ^z/^x and Bz/Sy. The ball rolling about the bowl may be 
visualized as describing the path traversed by the missile tip relative 
to the center of gravity. 
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Inserting a friction force proportional to rolling velocity and 
aerodynamic moments proportional to angular velocities results in 


X = K (x,y) - Lx 
ox 

y = K (x,y) - Ly 

oy 


for the ball, and 

a = Qm (a-,3) - lia 

3 = K 2 Cn (a, 3) - 

for the missile. In the case of a missile free to plunge, as well as to 
pitch, the coefficients contain several terms, but the differential equa- 
tion is in the same form as that presented. 

This, then, is the analogy between a ball rolling on a surface of 
given contour and a corresponding missile pitching and yawing with the 
above relations between surface slope and aerodynamic moments. 


EXAMPLES 


A case representing a missile with linear moment characteristics was 
tested to verify the analogue, and three nonlinear cases were tested, to 
demonstrate the application. All four cases were chosen axially symmetric 
to economize on model construction. The feature of axial symmetry sug- 
gests use of cylindrical coordinates to define the moment characteristics 
of the missile and the surface of the bowl. 

The moment tending to realign the missile axis with the flight path 
for the cases tested is given by 


Ccr = Ca + Da® 
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and Cp is related to Qm. and Cq. by 


Cin = cos e 
Cb = C 0 sin 0 


The bowl-surface slope is given by 


~ = 2Ar + 4Br3 
dr 


and the equation of the bowl surface is then 

z = Zq + Ar2 + Br'^ 

Converting the bowl-surface equation to Cartesian coordinates gives 

2 = Zq + A(x^y^) + B(x2+y^)^ 

The slopes of the bowl are then 

— = 2Ax + toxS + 4Bxy2 
dx 

— = -2Ay + 4By® + 4Byx^ 
dy 

Using the analogy between the bowl and the missile moment characteristics 
gives 

Cm = Ca + Da® + DoP® 

Cn = Cp + DP® + Dpa® 


These equations illustrate that Cm, for instance, is a function of both 
a and p, even for this simple nonlinear system. 
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The simplification attained "by use of cylindrical coordinates may 
apply to conventional electronic analogue computers as well as to the 
rolling-ball analogue. ^ 

The four cases were obtained by assigning various positive and nega- 
tive values to the coefficients C and D. The moment equations and the 
resulting bowl profiles are presented in figure 1. The terms "stable" 
and "unstable" in this figure refer to conditions near the bowl center. 
The cubic case has zero stability near the center of the bowl. 


CONSTRUCTION OF MODELS 


Several techniques of construction were tried. The most accurate 
bowl, also the least expensive, was cast in hard plaster over a die made 
of piaster. The first step in making a model was to make the die on a 
flat work surface by sweeping a steel ten^jlate about a vertical axis 
through the center of the die. While this die was hardening, a round 
frame was made using a piece of plywood for the base and a strip of sheet 
metal for the side. The hardened male die was greased lightly and the 
frame, with access holes cut in the base, was placed upside down over it. 
Plaster was then poured in through the access holes. The final contoured 
surface was then sanded and sprayed with flat, black lacquer. About eight 
man-hours of work were required per bowl for the axially symmetric cases 
built. These bowls were 2 feet in diameter. 


TESTING 


The bowls were placed directly below a motion-picture camera. Lag 
screws in the plywood base permitted accurate leveling of the bowls. The 
balls were all solid and were made of steel, brass, and aluminum and 
varied in size from 0.25 inch to 1.00 inch in diameter. No elaborate 
technique was required for starting the balls rolling, but it was found 

■^The equations of motion of the ball and missile may be converted to 
cylindrical coordinates. The equations for the ball are, assuming 
axial symmetry 

r _ re^ = K - Lf 
dr 

2re + re = - Lre 

The term 2re is the Coriolis force. The equivalent equations for 
the missile ara 


a 


ae^ = kQy - 


la 


269 + a’o = - laO 
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convenient to roll the balls down a chute rather than roll them by hand 
because good repeatability of initial conditions was desired. The motion 
pictures, which were teiken at either 8 - or I 6 - frames per second, were 
projected frame by frame onto a screen for measuranent. The camera fram- 
ing rate was used -as the time standard. Oil layers of various thick- 
nesses on the bowls were fo\md to provide good control of the dangoing. 
Using balls of various densities also provided scaae control over dan 5 )ing 
without oiling the surface. Actually, the damping was measured rather 
than predetermined. 


RESULTS 


The results obtained with the linear-case bowl are easily compared 
with analytic results and serve to check the validity of the analogue. 

In figure 2 are plotted the data for three runs corresponding to simple 
pitching motion with progressively heavier damping. This large variation 
in damping was attained by using a 1-inch-dlameter brass ball on the bowl 
with no oil for the light damping (fig. 2(a)) and a l/4-inch-dlameter 
aluminum ball with a thin coat of oil for the case of heaviest damping 
(fig. 2(c)). Superposed on each set of data is a damped sine wave fitted 
by the method of least squares (ref. l) . That the experimental data fit 
the theory quite well is clear. The undamped natural period, as predicted 
on the basis of bowl profile and the inertia characteristics of the solid 
balls used, is indicated in figure 2. 

Several plots of angle of attack versus time for pitching missiles 
with nonlinear moment c\irves, obtained with the bowls, are presented in 
figure 3« The curves for the cubic case (fig. 3(a)) were obtained at 
different times during the same test and Illustrate the extreme depend- 
ence of frequency on amplitude of oscillation, a characteristic attribu- 
table to the nonlinear moment curve. Ihe large -amplitude curve has the 
sharp peaks typical of cubic systems. The stable cubic case (fig. 3(b)) 
exhibits these same characteristics to a lesser extent. The unstable 
cubic case (fig. 3(c)) exhibits 6 points of Inflection per cycle instead 
of the 2 per cycle shown. by all the other cases. 

Two plots of motions obtained with the unstable cubic case are shown 
in figures 4 and 5 . The motion of figure 4 gradually degenerates and 
takes on a form much like that of figure 5* These two motions were not 
obtained on the same test run. 

In the course of experimental firings of a simple body of revolution 
in the Ames supersonic free-flight range, the pitch-yaw records shown in 
figure 6 were obtained. In the absence of large gyroscopic forces (the 
models were fired. from a smooth-bore gun) , it was difficult to explain 
this behavior. In attempting to do this, the analogy between the ball 
and the missile was discovered and the improved intuitive feel for the 
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problem made the explanation clear. The moment curve for this body has 
not been determined, but the motions observed indicate that the unstable 
cubic case is a fair representation. The data were obtained using bal- 
listic "yaw cards." The technique consists of deducing the model orien- 
tation at any instant from the shape of the hole it punches in a piece 
of paper. Therefore, the curves are not exactly defined, but the general 
nature of the motion is illustrated. 


LIMITATIONS 


The most serious limitation on the use of the analogue is that only 
pitching and yawing motion may be simulated. The missile which is being 
simulated must be assumed to be roll stabilized to zero rate of roll. 

The damping coefficient must be the same for pitch and yaw, which also 
limits the value of the analogue for airplane-like cases. 

It was noted diiring the test with the linear exan 5 >le that when the 
ball was executing elliptical motion, the axes of the ellipse precessed 
slowly. This distortion of the motion is believed to result from second- 
order momentum effects (arising from the fact that the bowl surface had 
finite slope) not considered in the derivation. The salient features of 
the motion are not masked by this effect. 


CONCLUDING KEMBKS 


The use of a ball rolling on a suitably shaped surface to represent 
the pitching and yawing motion of a missile has been described. Intuitive 
\mderstanding of some ccsirplicated motions is simplified throu^ the use 
of this device. The use of the analogue for axially symmetric cases sug- 
gests that the use of cylindrical coordinates for describing nonlinear 
moment characteristics may prove helpful in solving problems on conven- 
tional analogue computers. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Dec. 10, 1953 
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Moment curves and surface profiles for analogue 
surfaces tested. 
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Figure J.- Angle-of-attack histories of nonlinear cases, 
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(a) x^g A “ 0.464 



Angte of yaw, 

(b) Xcg /I ' 0.485 


Figure 6. - Experimental pitch-yaw histories for a body 
of revolution in free flight. Mach number • 45 . 
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